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We propose a theoretical model of quasi-one-dimensional superconductors, with attractive 
electron-electron interactions dominant in the singlet d-wave channel and sub-dominant in the p- 
wave channel. We discuss, in the mean field approximation, the effect of a magnetic field applied 
perpendicularly to the direction of the lowest conductivity. The lowest free energy phase corre- 
sponds to a singlet d-wave symmetry in low fields, but to a triplet symmetry in high fields. A first 
order singlet-triplet phase transition is expected at moderate applied fields of a few teslas. We pro- 
pose to ascribe the recent critical field and NMR experimental data, observed in superconducting 
(TMTSF) 2 C10 4 to such an effect. 

PACS numbers: 74.70.Tx, 74.25.Op, 74.20.Rp 
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The nature of superconductivity in the family of the 
quasi- ID organic superconductors ^MM (TMTSF) 2 X 
(X = PF§, CIO4,...) has been a long standing issue for 
the last three decades. There is still a debate whether it 
is a conventional superconductivity, with a completely 
gapped Fermi surface, or an unconventional one with 
points or lines of nodes on the gap. Many experiments 
have tried to address this question using different tech- 
niques such as NMR relaxation rate, thermal conductiv- 
ity, non magnetic impurity effect on T c , Knight Shift and 
upper critical field measurements. 

The temperature dependence of the proton spin lat- 
tice relaxation rate, measured by Takigawa et al. [J] in 
the CIO4 compounds, suggested the presence of lines of 
nodes in the gap on the Fermi surface. On the other 
hand, Belin et al. [5| measurements of the temperature 
dependence of the thermal conductivity argued for a fully 
gapped Fermi surface. More recently, Joo et al. [f| stud- 
ied the effect of non magnetic impurities on the super- 
conducting critical temperature T c 0], and showed un- 
equivocally that the order parameter changes sign on the 
Fermi surface and that consequently the gap has nodes 
where it is zero on the Fermi surface. This is in complete 
agreement with a model of unconventional superconduc- 
tivity in these organic materials. 

However, whether the pairing of the electrons in the 
superconducting state is in the singlet or the triplet sym- 
metry remains an unsolved question. Measurements of 
the upper critical field of superconductivity in these ma- 
terials [1, 0, [l(| showed a superconducting state surviving 
up to a field as high as 9 T in the case of the PFq com- 
pounds [a] , and at least 5 T in the case of the CIO4 com- 
pounds :9| . These fields exceed by far the Clogston-Pauli 
paramagnetic limit 11] for homogeneous singlet super- 
conductivity, which is estimated to be of the order of 
1.84 T c , with a T c around 1.1 K in these materials. This 
result is a strong indication that at high magnetic fields 
superconductivity could not be of a homogeneous sin- 



glet type. Many theories [H, 0J, 03 have been proposed 
to explain this behaviour by a triplet superconducting 
state, which is a non Pauli limited state. At low mag- 
netic fields, however, besides the observation by Andres 
et al. [15] of a diamagnetic signal when the field is per- 
pendicular to the chains indica ting a singlet state, mea- 
surements of the Knight-Shift [16[ showed a significant 
variation from the normal state behaviour of the spin 
susceptibility of Cooper pairs which decreased with de- 
creasing temperature in the two directions of the a— b' 
planes. This behaviour was interpreted as being the sig- 
nature of a singlet state at low fields. Therefore, it seems 
that there is a contradiction between a low magnetic field 
singlet behaviour and a high magnetic field non singlet 
one. 

Other authors, however [Tt], EH, pointed out that a sin- 
glet superconducting state may survive the Pauli param- 
agnetic limit by allowing for a spatially modulated order 
parameter, form ing an inhomogeneous singlet state, or 
FFLO state [l^, Ho, Hl[ , where the two electrons of the 
Cooper pair will have different wave vectors to compen- 
sate for the Zeeman splitting. Maki et al. was the first 
to consider FFLO in d-wave superconductors in CeCoIn^ 
and to consider the orbital effect [22I [23j . Nevertheless, 
a phase transition to an inhomogeneous FFLO state has 
not yet been identified experimentally with certainty in 
Bechgaard salts. 

Another important property of these quasi-lD materi- 
als, which must not be ignored to explain their behaviour 
in magnetic field, is that under an increasing magnetic 
field parallel to the b' direction the electron motion gets 
confined to the a— b planes of the crystalline lattice: when 
the magnetic field is large enough (H > 2T), the elec- 
tronic orbitals in each a— b plane are decoupled. There- 
fore the orbital pair breaking effect is almost suppressed 
and can be neglected in that case [U, Using this 

particular characteristic, we study a simple model which 
considers the possibility of the phase transition between 
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a singlet superconducting state at low magnetic fields, 
and a triplet one at higher fields, in order to explain the 
discrepancy between the low field singlet behaviour and 
the high field triplet one. 

The purpose of this article is to describe this transition 
and not to give the most accurate and most realistic val- 
ues of the critical fields H C 2 (T) , as we have neglected the 
orbital pair breaking effect. Nevertheless, it is true that 
this orbital effect is much weaker in Bechgaard salts than 
in isotropic metals, because of their lamellar character, 
as discussed below. 

Eventually, the transition between singlet and triplet 
superconductivity results from the small difference be- 
tween the corresponding values of the renormalized scat- 
tering coefficients g B and gt , as they have been calculated 
by a renormalisation group method[26l. |27|. This implies 
indeed a small free energy gap between both configura- 
tions. 

Theoretical model One of the most important features 
of the (TMTSF) 2 X compounds is the very weak elec- 
tron transfer in the c direction transverse to the a— b 
planes 0, [28| . The transfer integral in this direction is 
t c = 10 K, which is an order of magnitude less than 
that in the b direction, with tf, = 300 K, and two or- 
ders of magnitude less than that in the chains direction 
a, with t a — 3000 K. The particular structural feature 
of these compounds gives them a lamellar character, and 
hence reduces considerably the orbital destructive effect 
against superconductivity. Furthermore, when the mag- 
netic field is oriented parallel to the b' direction at low 
temperatures, the semi-classical electron motion can be 
shown to be oscillatory in the c direction [2^| with the 
amplitude 



Sz = 



4t e 



e VfH 

where Vf is the Fermi velocity, e the electron charge and 
H the magnetic field. 

Such a mechanism is somewhat analogous to the 
field-induced one-dimensionalisation leading to the 
F1SDW 0, ElO]. The electron system, in fact, under- 
goes a dimensional crossover from a low field 3D regime 
to a higher field effective 2D regime, where the elec- 
tron motion is confined to the a~b planes and the in- 
tcrlayer hopping is considerably reduced. This phe- 
nomenon has first been predicted by Lebed [24[ based 
on a quantum mechanical calculation taking into account 
the highly anisotropic shape of the Fermi surface. Magne- 



toresistance measurements along the c direction 
have confirmed these predictions. This confinement phe- 
nomenon reduces the screening currents in the c direc- 
tion and therefore cancels out their destructive effect on 
superconductivity. In our model, we take these consid- 
erations into account by assuming that the orbital pair 
breaking effect of the magnetic field can be ignored as 
a first approximation. Therefore, the only pair break- 
ing effect of the magnetic field that will be taken in our 
case is the Zeeman splitting of the electron energy. The 



inter-plane transfer, t c , is implicit, so that the mean field 
approximation is justified. However, it is weak enough 
to be ignored in the subsequent calculations. 

We consider a mean field Hamiltonian which allows 
for both singlet and triplet order parameters. At high 
magnetic fields the electron spin susceptibility is that 
of the normal state. This is why we assume that the 
field induced triplet state is the equal spin pairing state 
(ESP) [34], for it has no S z = component and therefore 
has the same spin susceptibility as the normal state. We 
will consider it throughout this paper. Furthermore we 
assume a unitary state with real order parameters. The 
singlet and triplet order parameters are assumed to have 
(i-wave and p-wave orbital symmetry respectively, see 
equations (5) and (6). The FFLO case, where the mixed 
singlet and triplet order parameters are spat ially modu- 
lated, was previously studied by Shimahara |35l . l36j . In 
our model we will only consider the homogeneous state. 
A state in which singlet and triplet symmetry coexist 
does not seem favourable in our model. 

The mean field Hamiltonian is described as follows 



H = H a +H s +H t 



H a = y £,r ct c r 



(1) 



(2) 
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H s = £ Mk){c_ %L% + 4 T c^} (3) 
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with the following definitions for the order parameters, 
where s and t subscripts stand for singlet and triplet 
respectively, 

A s (fc) = A s , n (fc) (5) 



A t (fc) = A t)TT (fe) = -A t)U (fe) 

= 2^T4(fc,fe')«c_p T c s , T }) 



(6) 



with 



V s (k,k') = V s °(2(k.k') 



- 1 



(7) 

V t (k,k') = V t °k.k' (8) 
In our model, we use the following assumptions : 

A s (fc) = A° S (2P X - 1) (9) 
A t (fc) = A° t k x and A t , u (fc) = (10) 
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V s (k, k) and Vt(k, k ) are respectively the d-wave and 
p-wave parts of the interaction potential, which is as- 
sumed to be attractive in both parts. The possibil- 
ity of coexistence of these two attractive parts of the 
interaction has been previously discussed by many au- 
thors (H, [13] . In the g-ology model @, HI] both interac- 
tions are included in the Hamiltonian and their relative 
relevance determines the symmetry of the superconduct- 
ing state. In our model, we do not discuss the physical 
origin of these interactions, whether they are mediated 
by phonons or by magnetic fluctuations, or any other 
mechanisms. We assume the coupling constants, V s and 
Vt, to be temperature and magnetic field independent. 
On the other hand, in the case of a very strong magnetic 
field, the coupling constant may indeed depend on the 



field strength. This case, which has been discussed by 
Kuroki et al. [39(, will not be discussed here. 

As a first approximation we will consider an isotropic 
dispersion relation. We think that the anisotropic aspect 
of the dispersion relation will not change qualitatively our 
results, as will be reported elsewhere. In the following, 
we will study the relative stability of a (i-singlet and a 
p-triplet superconducting state. We are interested, in 
particular, in the possibility of a phase transition between 
these two states in an increasing magnetic field. 

Self consistent equations It is not difficult to diago- 
nalizc the model Hamiltonian eq. (1). Here, we use the 
thermal Green's functions to solve it. The equations of 
motion for the relevant Green's functions yield the fol- 
lowing equations in matrix form: 



^s,k 

A,* 



A,s 



o 



(ii) 



where u> n = (2n + l)irT. Solving eq. (9) for the relevant A S (T, H) is given by 



Green's functions, ((etc* -,)) and ((cL c T - )), and per- 

forming the frequency sum we can write the self consis- 
tent equations for the order parameters of equations (5) 
and (6) as 



k 



{{&l P +{n e Hy-E c%l )thc- 

-{A 2 t:kl + {^Hf + E c:k ,)thc + } (12) 
_ V t (k,k')A. p 

- e ;/ - We* - Kp + ^ thc - 



c.k' 



(E cP +A 2 R ,-fi e H0 thc+} 



(13) 



where y, e is the electron magnetic field and we have de- 
fined 



E 2 t = £p + A% 

o,k k s.k 



(14) 



o,k s.k t.k 



VJk,k')A ,-, Q 
\,i = ~E 9B; /' fc {tanh(f(^- Me g)) 



IE ,-, 

w °' k 

+ tanh(|(Sg / + / x e ff))} 



(16) 



which is nothing but the BCS gap equation with the mag- 
netic field acting only on the spins of the electrons. When 
A s = in equation (13) we get another correct limiting 
case, that is, 



A a = -Ew')A a , 



t^mj^-^Hy + AK, ) 



(t«-»eHy + AK k> 



(17) 

This determines the field and temperature dependence 
of the p-triplet order parameter. 

By numerically solving equations (16) and (17) for A° 
and A° with the d and p symmetries given respectively 
by equations (9) and (10), we could calculate the critical 
field of superconductivity in these two limiting cases, us- 
ing the free energy difference with the normal state given 
by 



tanh(£ /££ + {^ e Hf±2E cP ) 

thc± = V (15) 

El + (/i e iJ) 2 ± c, k> 

It is straightforward to see that if we put A t = in 
equation (12) the field and temperature dependence of 



8F 







(18) 



Using equation (18) we could also calculate the critical 
field beyond which the <i-singlet state becomes energeti- 
cally less favourable that the p-triplet one. In fig. 1 are 
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Figure 1: Critical fields of superconductivity in the d-singlet 
(solid line) and p-triplet (dotted line) case, with a triplet criti- 
cal temperature of (a) 0.62 K and (b) 0.22 K. The dash-dotted 
line is the transition line between the two phases. The infinite 
slope of the triplet linte is unrealistic, and a consequence of 
the neglectance of orbital pair breaking effect. 



shown the phase diagrams of the pure d-singlet and p- 
triplet cases, with the transition line between these two 
phases. The strength of the triplet interaction, and ac- 
cordingly the triplet zero-field critical temperature, are 
unknown. Therefore we take the triplet critical tempera- 
ture as our free adjustable parameter. What should also 
be noted is that this transition to the triplet phase is a 
first order transition. However, accurate enough experi- 
mental measurements of thermodynamic quantities, such 
as specific heat or magnetic susceptibilities, have yet to 
be done to address the existence or the order of such a 
phase transition. 



Let us now study the possibility of a coexistence phase, 
where both the singlet and triplet order parameters are 
simultaneously non zero. When solving numerically the 
two coupled equations (12) and (13), we could not find 
any physically plausible solutions, which makes the coex- 
istence state very unlikely to happen. To confirm this re- 
sult we calculated a Ginzburg-Landau free energy devel- 
opment to the fourth order in A s and A t , see eq. (19). In 
the diagrammatic expansion of the free energy there are 
only three non-zero diagrams contributing to the A s 2 A t 2 
term. They are shown in fig. 2. All three diagrams have a 
positive sign in the T—H plane, which makes coefficient 
c in eq. (19) positive for all values of T and H, making 
the mixing term energetically unfavourable. 



SF C = SF + aA 2 



a' A 2 



6A 



b'At 



cA 2 A 2 (19) 



with 8F C being the free energy difference of the coexis- 
tence phase to the normal state. 

Discussion In conclusion, we have presented a theo- 
retical model in which we assume the simultaneous ex- 
istence of singlet and triplet pairing interactions. This 
model, treated in a mean field approximation, shows a 
competition between two homogeneous superconducting 
phases, a d-wave singlet phase, stable in low applied field 
and a p-wave triplet superconducting phase, stable un- 
der large applied field. We have shown the existence of a 
phase transition from the d-singlet phase to the p-triplet 
phase in a moderate field, of the order of two teslas. A 
phase in which coexist these two symmetries is not stable 
in our simple model. 

In this model, we therefore predict a very large critical 
field, much larger than the Clogston-Pauli limit. We are 
also able to understand the very peculiar behaviour of the 
NMR data recently observed in {TMTSF) 2 ClOi 0: a 
usual Meissner effect is expected in the singlet phase, 
while one should recover the normal phase susceptibility 
in the Equal Spin Pairing triplet phase, as experimentally 
observed. 

A qualitative property concerning the anisotropy of the 
critical field can be easily predicted. As the applied mag- 
netic field is rotated in the a—b plane, the field orien- 
tation corresponding to the highest critical field should 
depend on the field strength. In low fields, typically H 
smaller than about 2 teslas, the orbital effect of the field 
is too weak to induce the two-dimensional confinement 
observed at larger field [29]. Therefore, in this low field 
limit, we expect that the orientation of the highest crit- 
ical field, which, of course is determined by the orbital 
effect is along the a-direction, in which the electron trans- 
fer integral is the largest. On the contrary, in the large 
field limit, the field induced two-dimensional confinement 
strongly reduces the electron transfer in directions per- 
pendicular to b' . We therefore expect, as the field is in- 
creased above about 2 teslas, that the orientation of the 
maximum critical field is progressively rotated from the 
a-direction to the 6'-direction. 

In the large magnetic field case, the plane a—b are de- 
coupled and the orbital pair breaking effect is drastically 



Figure 2: Diagrammatic expansion of the A s 2 At 2 term in the free energy. 



reduced and can be neglected, as we have already dis- 
cussed. In the small magnetic field case, however, we 
expect that the remaining orbital effect will be similar 
for the singlet and for the triplet cases, so that the differ- 
ence between the triplet and the singlet free energy is not 
noticeably modified; therefore we hope that the orbital 
corrections will not destroy this transition. Moreover, we 
observe that experimental values of the critical field H c a 
are noticeably increased when the field is aligned along 
the b' direction; this effect is coherent with our assump- 
tion of a two-dimensional nature of the singlet/triplet 
transition. 

We have examined the case of a p-wave triplet state. 
We could have extended our study to an /-wave triplet 
state, which mi ght be favoured by two-dimensional fluc- 
tuations 0, 126ll27l |4C| . We believe however that similar 
results would have been obtained. 

Our two dimensional model, described in eq. (1-4), can 
be obtained by a three cut-off renormalisation group ap- 
proach [4l[ , in which the three cut-offs are the bandwidth 
Eo ~ 6000K, the cross-over temperature T co w 100K and 
the Zeeman energy fj, e H w 2K, and have different orders 
of magnitude. In the range where to < T co , our Fermi 
liquid analysis is valid. Since [i e H < T co , the Zeeman 
energy can also be discussed in a mean field approach, as 
we did here. 

Many authors [H, [H, [H, |42j have proposed to in- 
terpret the experimental data by the field induced sta- 



bilisation, not of a triplet phase, but of a singlet-FFLO 
phase. Both phases have somewhat similar properties : 
they appear in a first order transition ; they are very 
sensitive to non magnetic disorder, which might be the 
case of the high field phase ; they have a critical field 
which exceeds the Clogston-Pauli limit. However, there 
is a second paramagnetic limiting field beyond which the 
FFLO inhomogeneous singlet superconducting state will 
also be destroyed. This new paramagnetic limit was cal- 
culated by Lebed [3 in the quasi-lD (TMTSF) 2 X su- 
perconductors along the b' direction, the direction of the 
highest critical field, and was estimated to be given by 
Hp fflo = 0.6(t Q /ih) 1 / 2 i?p , where H p is the usual Pauli- 
limited critical field, t a and % are the electronic transfer 
integral along a and b. In the Bechgaard salts iTpFFLO 
should be of the order of 4 T. This limit was calculated 
without considering the orbital pair-breaking mechanism, 
which should reduce it further. It is not obvious to us 
that such a theoretical estimation of the FFLO critical 
field is large enough to account for the experimental val- 
ues. 

It seems to us that more experimental work is still 
necessary to discriminate between the various theoretical 
models proposed to interpret the experimental data 
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